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ABSTRACT 

Globular clusters (GCs) and many nuclear clusters (NCs) show evidence for hosting multiple generations of 
stellar populations. Younger stellar populations in NCs appear to reside in disk like structures, including the 
nuclear cluster in our own Galactic center as well as in M3 1 . Kinematic studies of the anomalous globular 
cluster cl) Centauri, thought to possibly be a former dwarf galaxy (or a galactic nucleus), show evidence for 
its hosting of a central, kinematically cold disk component. These observations suggest that formation of 
second (or multiple) generation stars may occur in flattened disk like structures. Here we use detailed N-body 
simulations to explore the possible evolution of such stellar disks, embedded in globular clusters. We follow 
the long term evolution of a disk like structure similar to that observed in co Centauri and study its properties. 
we find that a stellar-disk like origin for second generation stellar populations can leave behind kinematic 
signatures in properties of the clusters, including an anisotropic distribution, and lower velocity dispersions, 
and we discuss their implications in constraining the origin of second generations stars and their dynamical 
age. 



Subject headings: methods: numerical, 
galaxies: nuclei 



Galaxy: globular clusters: individual (w Centauri), Galaxy: disk. 



1. INTRODUCTION 

Nuclear star clusters (NCs), stellar systems similar to glob- 
ular clusters (GCs), have been found in ~ 75% of the late 
type spiral galaxies and dwarf elli pticals (Boker et al.l l2004l : 
IWalcher et al . 2005: Cot e, R et al. l l200 6: Boke r 2010). They 
contain several distinct stellar populations, characterized by 
different ages and chemical properties. Among these groups 
of stars, the youn gest ones often form centrally concen - 
trated stellar disks (IVan den BoschI 119941: ISeth et al.ll2006l) . 
Such stellar structures are also found in the central pc of the 
Galaxy and in M3L Their populations of massive young 
(6 + 2 Myr and ~ 200 Myr, respectively) stars a re ar- 
ranged in one (orrnore) distinct stellar disks (lT remaine 1ll995l: 
Genzel et al.l2000tlLevin & Beloborodovl2003i:iBender et al.l 



2005HLU et alJl2009l:lBartko et al.ll2010l:lGenzel et al.ll2010t) . 

Similarly, and in contrast with previous paradigm, GCs are 
not composed of a simple stellar population born at the same 
time from the same gas. Rather, stars in GCs show signifi- 
cant inhomogeneities in the abundance of light elements and 
their color magnitude diagrams show several isochrones; thus 
they consist of multiple stellar generations (MSGs) whose age 
spr eads is estimated to extend up to a few hundred Myrs (see 
e.g. lGratton et~ani200ll20T2h . 

The origin of these MSGs is still debated; in one of the 
models the low velocity gas ejected by primordial asymp- 
totic giant branch (AGB) stars is retained inside the clus- 
ter, concentrate s at its center and there fragments forming 
new stars (see iGratton et al.l 120121 and reference therein). 
In this scenario, the initial distribution of second genera- 
tion stars strongly depends on the initial profile of the col- 
lected gas (a somewhat similar model was suggested fo r 
the origin of M31 young stellar disk by IChang et al.ll2007l) . 
For instance, if the gaseous material has some angular mo- 
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mentum, it could produce a flattened structure instead the 
commonly assumed spherical config uration (Bekki 20111 
iMastrobuono-Battisti & Perets 1120131) . in analogy with what 
is observed in NCs. This would influence significantly the 
later dynamical evolution and mixing of older and younger 
stars, as we show in the following. 

One of first clusters where MSGs have been found is 
CO Centauri (NGC 5139, o) Cen). With mass estimates 
between 2 x 10^ and 5 x 10*' Mq (see lMevlanetal.1 119951: 
Ivan de Ven et al.l l2006h . co Cen is the most luminous and 
massive Galactic GC . It is characterized b y a flattened 
structure (iGeyer et al.l 119931 Ivan de Ven eTa l. 2006 | and 
by a wide range of me tallicities (e.g. iButler et al.l 119781 : 
iNorris & Da Costa I [T995I) . Its MSGs are characterized b^ 
different spatial d istribu tio ns (e.g. [F reeman & Rodgers 



T97I iLee et al. 



Hilker & Richtler 



iBedin et al. 



1 1 999t iHughes & Wallerstein , ,200' 



20041: IStanford et al 



20001: jPancino et al.l 120001: iFerraro et aL 



2OOI iSolUmaetai 



Dynamically, o) Cen's orbit is peculiar, being 
tightly bound to the Galactic poten tial, retrograde and 
highly eccentric (iDinescu et al.l 119991) . Moreover, it is 
one of the favorite candidates for hosting an intermediate- 
mass black hole, whose mass is estimated to be in the 
range 13,000 < M. < 50,000 Mp ( Novolaetal. 200: 



van der Marel & Anderson 1120101 iMiocchi Il20ia IJalah et al 



2OI2I) . Given these features a) Cen is thought to be the 
remnant nuclei of a satelUte dwarf galaxy that merged 



with the MiUcy Way (iFreemanI IT993t Dinescu et al.l 
Hughes & Wallerstein I 120001: iBeldd & Freeman 



Bokerll2008l) ." 



1999 



2003 



Recentlv Ivan de Ven et akl (120061) used an extension of the 
ISchwarzschild I(ll979l) method to build axisymmetric dynami- 
cal models for w Cen, using anisotropic velocity distributions. 
Their best fitting model shows a highly isotropic inner region, 
which becomes increasingly anisotropic outside of 10 arcmin. 
The nucleus of w Cen contains a kinematically cold disk com- 



Mastrobuono-Battisti & Perets 



TABLE 1 
Initial conditions for the simulation runs. 



Label 


N 


m..i 


m,i 


Mdisk 


'V 


potential 






Mg 


Mo 


Mo 


pc 




Sl,2 


4000 


25 





10^ 


4.5 


isolated 


S3 


10000 


25 


12.5 


1.5 X 10^^ 


4.5 


isolated 


S4 


8000 


12.5 





105 


4.5 


isolated 


S5 


4000 


25 





105 


4.5 


tides 



Note. — SI and S2 are two different realizations of the 
same disk; S3 is a simulation of the disk with two mass 
classes but with the same relaxation time; we run also a 
case with a disk as massive as in Sl,2 but with a bigger 
number of stars (S4). Finally S5 is a simulation where 
the cluster orbits a point mass galaxy. 

ponent, with a radius between 1 and 3 arcmin, an average flat- 
tening of ~ 0.60, and contributing about ~ 4% to the total 
mass of the cluster 

The presence of this disk could be closely related to the 
existe nce of MSGs in cj Cen and inside other GCs (iBekkil 
1201 Ih . In particular, second generation stars may form in a 
disk, and leave behind observable kinematic signatures of 
their former configuration even after long term dynamical 
evolution. These signatures may allow us to infer information 
about their initial distribution, and provide us with constraints 
on their origin. 

In this paper we study the long term evolution of w Cen-like 
central disk by means of N-body simulations and explore its 
structure and potentially observable signatures. This test case 
is used as a pilot to reveal the possible kinematic signatures 
of second-generation stellar disks and point to possible clues 
for the initial presence of a second generation disk in GCs and 
NCs. 

The paper is organized as follows. In Section|2]we describe 
our choice of the initial conditions, the simulations we run and 
the code that we use. Then, in Section[3] we show the results 
obtained. These results and their implications are discussed 
and summarized in Section|4] 

2. THE N-BODY SIMULATION 

2. 1 . The simulation 

We studied the dynamical evolution of an o) Cen-like cen- 
tral disk structure by means of N-body s imulations. Thes e 
simulations have been done using 0Grape (iHarfst et al.l2006l) . 
a high precision, collisional, direct summation code originally 
configu red to mn on com puter clusters accelerated by GRAPE 
boards dMaldno & Taiii | [T998). This code has been recently 
adapted to run on GPU equipped machines by mean of Sap- 
poro, a CUDA library that calcu lates force on the GPU s mim- 
icking the behavior of GPAPEs ( iGaburov et al. Il2009h . 
The code implements a fourth-order Hermite integrator with 
a predictor-corrector scheme and hierarchical time stepping. 
The time-step parameter 77 and the smoothing length e are the 
parameters that controls the accuracy of the code. In our sim- 
ulations, we set rj - 0.02 and e = 5 x 10"~ pc. The code has 
been run on the cluster "Tamnun" at the Technion. 

2.2. Initial conditions 

The initial conditions for our simulations were set based 
on the available observational data. Since the relaxation time 
of w Cen is longer than a Hubble time ( .Harris, ,20 10,) , mass 
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Fig. 1 . — Projected isodensity contours plot of the disk on the X-Y plane, at 
the beginning of the simulation, and going downward, after 2, 5 and 12 Gyr. 



segregation is not expected to have had significant eff'ects on 
the system. Its structure is therefore modeled as a single- 
mass, spherical and isotropic .King , (,196^) cluster. To build 
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Fig. 2. — Cumulative distribution of inclination of stars in the disk at differ- 
ent, equally spaced, times. Going upwards the curves correspond to 0, 1.5, 3, 
4.5, 6, 7.5, 9, 10.5 and 12 Gyr 



our model we used the parameters that best fit the cluster ob- 
served brightness profile: we set Vc - 4.6 pc as the core radius. 
Wo = 5.5 ( iMeylan II 1987b and 2 5 x 1 0^ Mq as the total mass 
of the GC (see lvan de Ven et"ani2006l) . The N-body represen- 
tation of the spherical component consists of N - 100,000 
particles. 

In the main simulation we run (Sl,2), the disk mass is 
10^ M0, and, as the observed one, has a radial extension of 
4.5 pc and a flattening similar to that observed. For its N- 
body representation we used Nd - 4000 single mass particles 
orbiting counter clockwise in the X-Y plane. 

Since we are interested in the evolution of the very cen- 
tral region of to Cen we neglected the currently observed 
anisotropy of the cluster 



In addition, at least in the central region, the influence of 
the tidal field of the Galaxy is negligible as compared to the 
mutual interaction between particles, thus we simulated the 
cluster as isolated. Nevertheless, in order to exclude spurious 
effects caused by these assumptions we run another set of test 
simulations with different initial parameters for the disk (see 
Table [U. In particular, we put the cluster on a circular or- 
bit in a point mass potential as massive as the Galaxy within 
CO Cen's current position. Since the disk evolves in the same 
way as in the case of the isolated cluster, we can conclude that 
the presence of the galactic tide doesn't introduce any signif- 
icant effects. We also run a second realization of the same 
disk+cluster system and, in order to exclude mass segregation 
effects, we simulated a disk with two different mass classes 
(12.5 and 25 Mq) with the same relaxation time of the origi- 
nal model. In both cases the evolution of the system doesn't 
differ significantly from the first case. Finally, we tested the 
dependence on the chosen number of particles by simulating 
a disk with 2N particles (keeping the same total masses of the 
disks and cluster, but normalizing for the different relaxation 
time). Again, we find no significant changes compared with 
our main simulation. For simplicity, and given the insignif- 
icant differences between the various models we tested, we 
focus here on our main simulation (SI), and show only its re- 
sults, keeping in mind that they also well represent the results 
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Fig. 3. — The angular momentum of the disk. Upper panel: The evolution of 
the three components of the angulai' momentum of the disk particles. The L, 
Ly components lay on the initial plan of the disk, while L- is initially perpen- 
dicular to the disk. The disk loses angular momentum that is redistributed to 
the stai's in the cluster. Bottom panel: The dispersion of the three components 
of the angular momentum. 



of the other models (see TablefTli. 

3. RESULTS 

If second generation stars form in a disk, it is possible that 
their initial distribution leaves observable signatures in the 
parent dense cluster Thus, we simulated the evolution of the 
CO Cen-like central disk for a period of time that, re-scaled to 
the mass of the particles, is comparable with the typical 12 
Gyrs age of Galactic GC. We then analyzed the evolution of 
the angular momentum of the stars in the disk as well as their 
orbital parameters in order to find relevant kinematic proper- 
ties that might be detected in real clusters. 

Figure [1] shows the isodensity contours of the disk at the 
beginning of the simulation and after 2, 5 and 12 Gyr The 
disk inflates with time and becomes more spherical but its 
radial size is only slightly increased. The disk stars remain 
spatially limited inside ~ 15 pc, while the tidal radius of the 
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Fig. 4. — Top left panel: the velocity dispersion of the disk at four different evolutionary times, 0, 2, 5 and 12 Gyr; V; is the velocity component perpendicular 
to the initial plane of the disk. Top right panel: The final radial (blue thick and black thin soHd Hnes) and the tangential (blue thick and black thin dashed lines) 
velocity dispersion of the disk and of the cluster as a function of the projected radius. Bottom left panel: The final anisotropy parameter of the disk (solid line) 
compared with that of the cluster (dashed line), as a function of the projected radius. Bottom right panel: Relaxation time of the central region of the initial 
simulated system as a function of the radius. 



cluster is ~ 80 pc. 

The distribution of the average inclination of the disk stars 
becomes broader with time. After 12 Gyrs it is much more 
uniform (see the cumulative distribution evolution in Figure 
|2]), i.e. the disk becomes more isotropic with time, in 
agreement with the results in Figure [1] 

Figure |3] illustrates the evolution of the three components of 
the disk angular momentum (AM, upper panel), showing a 
consistent decrease with time, and the respective dispersions 
(bottom panel) that increase with time. At the end of the sim- 
ulation the average AM is still significant; even after the long 
term evolution the disk stars still do not show an isotropic 
distribution. Moreover, as the total AM of the system must 
be conserved, the cluster structure fl attens as it exchanges 
AM with the disk stars. Interestingly, Ivan den Berghl (12008b 
found a relation between the morphology of the horizontal 
branch and the flattening of Galactic GCs. It is possible that 



our result could be connected to this result, since the shape 
of the horizontal branch is known to be strongly related to 
the evolutionary history of the GC and therefore with the 
presence of MSGs. 

While the disk inflates, its velocity dispersion grows with 
time (see top left panel of Figure |4]i. Despite that, at the end 
of the simulation, the velocity dispersion of the disk stars 
is still smaller than that of the cluster stars. Both the ra- 
dial and tangential velocity dispersion of the disk and that of 
the cluster differ, with the former considerably smaller than 
the latter (see the top right panel of Figure |5}. Moreover, 
the disk remains significantly radially anisotropic in velocity 
space. Its anisotropy parameter, at the end of the simulation, 
/3 - 1 - o-j/{2cr;.) is indeed positive at any radius, and has an 
increasing trend with the distance from the cluster center (see 
bottom left panel of Figure |4|i. 
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Thus, even after 12 Gyrs the flat and the spherical compo- 
nents are not yet mixed and can still be spatially and kinemat- 
ically distinguished. This suggests that signatures of any pu- 
tative initial distribution of MSGs formed in a disk-like struc- 
ture might still be observable today. This is true at least in 
clusters (or NCs) with long relaxation times as <jj Cen, as well 
as in CO Cen itself, where many stellar generations have been 
found. Finally, as a result of the evolution of the disk, the 
whole system becomes tangentially anisotropic with time in 
the region outside 10 pc, in agreement with what found by 
Ivan de Ven et al, (,2006.) . 



4. DISCUSSION AND SUMMARY 

GCs are made of several generation of stars whose origin is 
still unclear, o) Cen is the first cluster where this feature has 
been observed. Among its peculiarities, this cluster shows a 
central disk that could be the remnant of the an initial disk 
where second generation stars formed. Indeed iBekki I (1201 lb 
showed that the gas ejected by AGB stars can be retained in- 
side the parent massive star cluster and can fragment into a 
nested flattened and rotationally supported system of young 
stars. If so, the flattened structure evolves and may eventually 
leave signatures of its former presence. 

In this work we employed N-body simulations to explore 
the long-term evolution of such stellar disks embedded in 
dense stellar cluster, where for the initial conditions we made 
use of the observed disk like structure in w Cen . 

We verified that such a disk becomes more isotropic with 
time. However, even after a time comparable with the age of 
Galactic GCs, the cluster stars and the second generation disk 
stars are not yet completely mixed. The two populations have 
different spatial distributions, with the second generation stars 
still confined only to the central region of the cluster, and do 
not attain a completely relaxed spherical shape at the end of 
the simulations. As the second generation disk stars isotropize 
they exchange angular momentum with the first generation 
cluster stars, which consequently leads to the slight flattening 
of the host cluster. It is worth noting that a more massive disk, 
as e.g. shown to form in simulations (.Bekki,.201ll) is likely 
to introduce an even larger flattening. 



We also find that the kinematic properties of the evolved 
disk stars change with radial distance. The disk preserves a 
slight (y6 ~ 0.05) radial anisotropy within 4 pc; outside this ra- 
dius p grows significantly. At the same time the whole system 
comes to be tangentially anisotropic out of 10 pc, in agree- 
ment with observations. In addition, the second generation 
stars are also characterized by a significantly smaller velocity 
dispersion than the first generation stars in the central region, 
which decreases even further with the distance from the clus- 
ter center 

In cases where relaxation processes are less efficient, initial 
conditions are better preserved, as the cluster kinematics still 
hold their dynamical "memory". This explains why the kine- 
matics of the disk stars show stronger signatures of their orig- 
inal disk structure as a function of increasing radius, at which 
relaxation times ar e longer (see bottom rig ht panel of Figure 
ID and Figure 21 in Ivan de Ven et all 120061) : indeed, a similar 
distance dependent effect was suggested to constrain the ori- 
gin of young stars in the G alactic center (iPerets & Gualandrisl 
l2010tlMadigan et al]l2013l) . A clear implication/prediction of 
this issue is that second generation stellar populations in GCs 
and NCs with longer relaxation times, and/or further out from 
the center, are likely to show stronger anisotropics in their 
kinematics. Such regions/clusters would serve as the best test 
sites for identifying fossilizes signature of a disk/flattened pri- 
mordial structure. 

Finally, the observed signatures of such disks may allow us, 
to some extent, to determine the dynamical age of the respec- 
tive stellar population. Since longer evolved disk stars show 
weaker kinematic signatures, the level of anisotropy provides 
a dynamical age stamp for this stellar population (see Figures 
[T] |2] and [Hi, and might even be translated into the real age 
of the stellar population, after accounting for the host cluster 
relaxation time. 

Follow up studies exploring larger primordial disk struc- 
tures than those studied here, beginning with the initial condi- 
tions j ust after second generation star forma tion (e.g. iBekki I 
120111: iMastrobuono-Battisti & Perets Il2013h . will aUow us to 
further explore the generic effects of the long term evolution 
of second generation stellar disks in a wider range of condi- 
tions. 
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